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The bandstructure of the zinc- blende phase of A1N, GaN, InN is calculated employing the exact- 
exchange (EXX) Kohn-Sham density-functional theory and a pseudopotential plane-wave approach. 
The cation semicore d electrons are treated both as valence and as core states. The EXX bandgaps of 
A1N and GaN (obtained with the Ga 3d electrons included as core states) are in excellent agreement 
with previous EXX results, GW calculations and experiment. Inclusion of the semicore d electrons 
as valence states leads to a large reduction in the EXX bandgaps of GaN and InN. Contrary to 
common belief, the removal of the self-interaction, by the EXX approach, does not account for 
the large disagreement for the position of the semicore d electrons between the LDA results and 
experiment. 
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I. INTRODUCTION 

The modern computational methods applied to con- 
densed matter systems are mainly based on Kohn-Sham 
PJ formalism of density-functional theory (KS-DFT). 
When used in conjunction with both the local-density 
(LDA) and generalized gradient approximations (GGA) 
for the exchange-correlation (XC) potential, the KS-DFT 
approach provides a very efficient and successful tool to 
compute the total energy and other related ground state 
properties. However, calculated KS bandstructures of 
solids display several limitations of this approach. The 
most serious shortcoming is the so-called bandgap prob- 
lem: the LDA bandgaps of semiconductors and insulators 
are between about 30 to more than 100% smaller than 
the corresponding experimental values. Another prob- 
lem is the underestimation of the binding energies of the 
semicore electrons. For the group-Ill nitrides and II- VI 
compounds, currently intensely studied materials due to 
their potential technological applications, the LDA re- 
sults for the position of the d bands are higher than ex- 
periment by 2 to 4 eV. The GGA electronic structures 
of these materials differ marginally from those of LDA. 
It should be noted that these deficiencies reflect on the 
one hand a fundamental issue: even the (unknown) ex- 
act single-particle KS eigenvalues will not directly corre- 
spond to excitation energies of the interacting electron 
system, except in special cases 0,0- As to the bandgap, 
it is well known [4| that there is a discontinuity in the 
exact XC potential (Axe) when an electron is added to 
the system. In general the magnitude of Axe is however 
unknown. The second, separate, issue is the effect of the 
(inevitable) approximation for the XC functional. The 
LDA errors for bandgaps as well as for binding energies 



of semicore states have been attributed [j| to a spurious 
self-interaction in LDA and GGA calculations 0, |(J . 

The best available technique for calculating the elec- 
tronic structure of solids is the quasiparticle theory, 
within Hcdin's GW approximation 7]. The non- 
sclfconsistent GW approach based on LDA results (LDA- 
GW) is found to reproduce very well the experimental 
bandgaps of standard sp semiconductors Q . For group- 
Ill nitrides and II- VI compounds, the pseudopotential 
HE3 and all-electron [n) LDA-GW methods are found 
to give good bandgaps (few tenths of an eV smaller than 
experiment), provided that the entire semicore shell is 
considered as valence in the pseudopotential approach. 
The positions of the d bands are markedly improved 
within LDA-GW USEI, although they remain system- 
atically higher than experiment by about 1 eV. 

A very interesting recent development is the applica- 
tion of the exact-exchange (EXX) approach |12f . within 
the KS-DFT approach, to calculate the electronic struc- 
ture of solids. In this approach the local - KS - ex- 
change potential is treated exactly, i.e., the spurious self- 
interaction [|| of the electron with its own charge density 
is eliminated. In the exchange- only KS-DFT formalism, 
the EXX potential is equivalent to the variationally 
best local approximation (i.e., the optimized effective po- 
tential) to the nonlocal Hartree-Fock operator 0, ^| . 
The EXX approach was first applied to crystalline solids 
by introducing a spherical shape approximation of the 
potential by Kotani and co-worker |16l 1 171 ll 8t [Tjj . This 
approximate EXX approach is found to highly improve 
the calculated bandgaps of several semiconductors 0] 
and the electronic structure of transition metals |19| , with 
respect to the LDA results. More recently, a pseudopo- 
tential EXX approach has been introduced |20i l2lj . The 
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pseudopotential approach is free from the above shape 
approximation of the potential and the convergence with 
respect to the computational parameters can be easily 
controlled. The pseudopotential EXX approach is found 
to give p]| bandgap energies in good agreement with ex- 
periment and GW results 0, H3 , for a large variety of 
semiconductors studied - except for diamond. On the 
other hand, for noble-gas solids 23] the pseudopoten- 
tial EXX approach gives bandgaps which are larger than 
those of LDA, but they are still much smaller than the 
experimental values. Thus, in addition to their own fun- 
damental and practical importance, the EXX results pro- 
vide a very good starting point for the non-sclfconsistent 
GW (EXX-GW) calculations Q. 

The purpose of this work is of threefold: (i) To perform 
pseudopotential EXX calculations of the electronic band- 
structure of the meta-stable zinc-blende (ZB) phase 
of GaN, A1N, InN. In these calculations the Ga 3c? and 
In 4d electrons are treated as valence states. The elec- 
tronic structure parameters of the ground state wurtzite 
(W) phase of the considered nitrides can be inferred from 
those of the corresponding ZB structure [26l |2*7| . (ii) To 
systematically analyze the effects of the semicore states 
on the bandstructure. To do that, an additional set of 
EXX bandstructures for GaN and InN are calculated, 
with the semicore states included in the core, (hi) To 
analyze and separate the changes that are due to the use 
of EXX pseudopotentials (PPs) from those due to the 
EXX treatment of the valence electrons. This is done 
by performing EXX bandstructure calculations employ- 
ing LDA-PPs, and LDA calculations using both LDA- 
and EXX-PPs. 

It is worth noting that the bandgap of WTnN is still 
a matter of debate. The bandgap of the recently grown 
single crystalline hexagonal InN epilayers is found |28j to 
be of 0.7-1.0 eV, which is much smaller than the previous 
and widely accepted |2{| experimental value of 1.9 eV. 
However, the smaller bandgap is far from being well es- 
tablish ed J30L l3l| . A very recent experimental study sug- 
gested 30] that this small bandgap is a misinterpretation 
of the Mie resonances due to scattering or absorption of 
light in InN samples containing clusters of metallic In. 
Moreover, a tentative bandgap of 1.4 eV has been given 
in Ref. y(J. Therefore, accurate theoretical studies of the 
bandgap of InN are in order. 

The paper is organized as follows. Sec. II contains 
a description of the computational method. Our results 
are reported and discussed in Sec. III. A summary of our 
main results and conclusions is given in Sec. IV. 



II. COMPUTATIONAL METHOD 

To perform EXX band structure calculations, we im- 
plemented the pseudopotential EXX formalism |2lJ in 
the SPHIngX code [33], which employs plane wave ba- 
sis sets and ab initio pseudopotentials. As noted above, 
EXX bandstructure calculations are performed for only 



the ZB phase of the considered nitrides. This is done at 
the experimental lattice constants (4.37, 4.50 and 4.98 A 
for A1N, GaN and InN, respectively [33| ). We use the ex- 
perimental rather than the theoretical lattice constants 
(which can be obtained by minimizing the EXX or LDA 
total energy) to allow for a direct comparison with the ex- 
perimental bandgaps and to avoid shifts in the bandgaps 
due to deformation potentials. 



Table I: Pseudopotential parameters. The local part 
of the ionic pseudopotential is chosen to be one of its 
components, referred to by its orbital angular momentum 
(I). For the Ga +3 and In +3 LDA-PPs, reference energies 
for the Ga Ad and In 5d states of 20 and 25 eV are used, 
respectively. 
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Norm-conserving scalar-relativistic EXX- and LDA- 
PPs of the involved elements are generated by using the 
Troullier-Martins optimization method [9. The EXX- 
PP's are constructed as described in Ref. [35|. Then, the 
generated pseudopotentials are transformed to the sep- 
arable Kleinman-Bylander form [3q|. The ground state 
electronic configuration is assumed for all the atoms con- 
sidered. Following Ref. l26l only s and p components are 
included for N. The pseudopotential parameters used are 
listed in Table I. The transferability of the pseudopoten- 
tials has been carefully tested, and they are found to be 
free from ghost states |32| • The construction of the LDA- 
PPs and the testing of both the EXX- and LDA-PPs are 
performed by using the FHI98PP code [38| . 

In the LDA calculations we used the Ceperley- 
Alder |39] exchange-correlation data as parameterized 
by Perdew and Zunger Q. In the present EXX calcu- 
lations, the exchange energy and potential are treated 
exactly and the above LDA functional is used for those 
of the correlation. Hereafter, unless otherwise specified, 
this combination (exact-exchange plus LDA correlation) 
will be referred to simply as EXX. Brillouin zone integra- 
tion is performed on a regular 4x4x4 Monkhorst-Pack 
|40j mesh. The KS wavefunctions are expended in terms 
of plane waves up to a cutoff energy of 60, 65 and 70 
Ryd for A1N, InN, and GaN, respectively. To calculate 
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the independent particle polarizability, Xo (see Ref. l2lT) . 
plane waves up to an energy cutoff of 35, 45 and 55 Ryd 
are included for A1N, InN, and GaN, respectively. The 
energy cutoffs used for the EXX calculations for GaN 
and InN, with the semicore d electrons included in the 
core, are exactly as those of A1N. The k-point sampling 
and the energy cutoffs are carefully tested, and are found 
to give an excellent convergence. For example, reducing 
the Xo energy cutoff of GaN from 55 to 45 Ryd, keeping 
the other cutoff fixed at 70 Ryd, changes the calculated 
eigenvalues by less than 0.01 eV. The calculation of the 
exact exchange potential involves also the unoccupied KS 
states, through xo and the Green's functions of the or- 
bital shifts. To ensure convergence, we fully diagonalizc 
the Hamiltonian in the plane wave representation and 
include all conduction states. 



III. RESULTS AND DISCUSSION 

In this section we report and discuss our EXX and 
LDA electronic structures of the ZB form of A1N, GaN 
and InN. We will also briefly address the modifications 
of their bandstructure when going from the ZB to W 
phases. The most prominent bandstructure parameters 
of the ZB phase of the three nitrides considered, obtained 
with the semicore d electrons treated as part of the frozen 
core, are listed in Table II. Whereas, those of GaN and 
InN, obtained with the semicore d electrons treated as 
valence states, are listed in Table III. The available ex- 
perimental results for GaN are also shown in Table III. 
As an example, we show in Fig. 1 the LDA and EXX 
bandstructures of ZB-GaN. 



A. Comparison between the electronic structures 
of ZB and W phases 

We focus here mainly on the differences in the bandgap 
and the position of the shallow d bands (Ed) by going 
from the ZB to W phases, studied at the LDA level. For 
the W phase, we have also considered the experimental 
structural parameters j3^|. The obtained results can be 
summarized as follows, (i) In agreement with previous 
LDA calculations [2(1 |2tJ , the direct bandgap at the T- 
point (Eg) of the W phase of GaN and InN are larger 
than that of the corresponding ZB phase by about 0.2 
eV. (ii) The ZB-A1N is an indirect bandgap semicon- 
ductor with conduction bands minima at the X-points, 
in agreement with experiment and previous calculations 
[M EH H3- (w) The values of Eg are almost identical in 
both ZB and W phases of A1N, in agreement with a recent 
LDA study 26] . (iv) The value of Ed remains almost the 
same by going from the ZB to W phases, for both GaN 
and InN. For instance, the LDA results for Ed of GaN 
(obtained using EXX-PPs), with respect to the valence 
bands maximum, are of -16.84 eV for the ZB phase and - 
16.87 eV for the W phase. In the following we will mainly 



concentrate on the EXX and LDA electronic structure of 
only the ZB phase of the considered nitrides. For this 
reason, the reference to the ZB phase will be suppressed, 
and that to the W form will be explicity stated. 



B. Bandgaps and the effects of the semicore d 
electrons 

We will start our discussion with A1N where no semi- 
core states are present. Table II shows that the EXX 
results for the fundamental indirect bandgap (Eg ) are 
in very good agreement with the experimental value 
(5.11 eV, see Ref. 0), and with the GW results @. For 
Eg, experimental results are only available for the W 
phase, and the most recent value is of 6.033 eV (Ref. |43|) . 
Since the values of Eg are almost identical in the ZB and 
W structures of A1N, one can conclude that the EXX 
and GW results for Eg are in good agreement with each 
other, and with the above experimental value. 

We turn now to the electronic structure of GaN, ob- 
tained without the semicore d states and using EXX-PPs 
(Fig. la). This figure shows that the energetic position 
and dispersion of the valence bands are weakly affected 
when going from the LDA to EXX. The most noticeable 
difference is in the upper valence bandwidth. The con- 
duction bands, however, are significantly shifted upward 
in energy in the EXX calculations, with respect to the va- 
lence bands maximum. This upward shift increases the 
bandgap from 1.97 eV (LDA calculation with EXX-PPs, 
see Table II) to 3.52 eV. This EXX bandgap is in excel- 
lent agreement with previous EXX calculations HHH, 
LDA-GW calculations 10] (with Ga 3d electrons were 
included in the frozen core) and experimental data (3.2 

- 3.3 eV, see Table. III). 

To analyze the effects of the semicore electrons, we 
show in Fig. lb the EXX and LDA bandstructures of 
GaN, obtained by treating these electrons as valence 
states. In this case a rather different behavior is found: 
Compared to the EXX calculation without semicore 
states, the bandgap closes by 0.64 eV — leading to a 
bandgap of 2.88 eV. This reduction can be easily under- 
stood as a consequence of the symmetry allowed anion p 

- cation d repulsion |44| , which pushes the upper valence 
bands states up in energy, leading to a smaller bandgap. 
This EXX bandgap is in remarkable agreement with the 
pseudopotential LDA-GW result (2.88 eV) of Rohlfing et 
al. [ljj, obtained with the entire Ga semicore shell con- 
sidered as valence. Furthermore, these pseudopotential 
EXX and LDA-GW results agree well with the mixed 
basis all-electron LDA-GW result (3.03 eV) for W-GaN 
[ll|. after subtracting the difference in E^ between the 
W and ZB phases of GaN (~ 0.2 eV, see above). This in- 
dicates that the inclusion of the entire semicore shell as 
valence is not necessary in the EXX and EXX-GW (see 
Ref. 124)) calculations. This is largely due to the linearity 
of the EXX potential j3f| and the improved description 
of the semicore atomic s and p states of Ga by the EXX 
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approach. The binding energies of these states are ap- 
preciably increased by the removal of self-interaction, as 
shown in Fig. 2. Compared to experiment, both EXX 
and LDA-GW approaches underestimate the bandgap of 
GaN (by about 0.4 eV). However, it has already been 
shown that a combination of these two methods (i.e., the 
EXX-GW approach) gives [24[ bandgaps that are in very 
good agreement with experiment for GaN and some II- VI 
compounds. 

Concerning InN, the experimental bandgap (available 
only for W-InN) is still controversial, as noted above. 
Our results listed in Tables II and III for InN show a 
behavior which is qualitatively very similar to that of 
GaN: Inclusion of the semicore (In Ad) states in the core 
gives an EXX value for Eg of ~ 1.5 eV, while consider- 
ing them as valence lowers significantly this bandgap to 
~ 0.8 eV. Assuming a similar error as that of the bandgap 
of GaN (of about 0.4 eV) and taking into consideration 
the difference in the bandgap of InN when going from 
ZB to W phases (see Sec. III. A), a direct fundamental 
bandgap for W-InN of ~ 1.4 eV is estimated. This re- 
sult is smaller than the previously accepted experimental 
value of ~ 1.9 eV, and larger than the very recent results 
(0.7-1.0 eV, Ref. H^). Incidentally, it agrees well with 
the (tentative) experimental value of 1.4 eV suggested by 
Shubina et al. 30] . However, it should be noted that the 
EXX-GW approach does not always give bandgaps that 
are larger than those of EXX (Ref. @), and, hence, the 
present estimate should be also treated as tentative. 

It is interesting to note that the EXX bandgap of InN 
when including the semicore (In Ad) states is significantly 
larger than in all-electron LDA-GW calculations [TlLl47| . 
which give a bandgap of ~ eV, for the W phase. This 
is in contrast to the case of GaN where both approaches 
give almost identical results. The reason appears ji]} to 
be related to the fact that the dielectric function (which 
enters the GW) is calculated using LDA, which gives 
a semimetallic ground-state (negative bandgap, see Ta- 
ble III) for InN. Thus, the calculated dielectric function is 
rather different from the real one which should be that of 
a semiconductor. It is well known that Ge has also a neg- 
ative LDA bandgap (about -0.1 eV). However, since this 
bandgap is smaller than that of W-InN and the overlap 
between the conduction and valence bands is localized 
in a very small region around the V point, this does not 
seem to affect the application of the LDA-GW method to 
Ge (see for example Ref. @) ■ To address this issue in the 
case of InN two indirect LDA-GW approaches have been 
adopted in the literature: (i) Bechstedt and Furthmiiller 
[27| have performed pseudopotential LDA-GW calcula- 
tions with the semicore d electrons included in the core. 
The reduction in the bandgap due to the inclusion of 
these semicore electrons as valence is taken into account 
a posteriori: assuming the same value as found by LDA 
calculations, (ii) In the work of Usuda et al. [43, the 
bandgap has been opened by applying hydrostatic pres- 
sure, i.e., by reducing the InN lattice constant. The ac- 
curacy of these indirect LDA-GW approaches is difficult 



to assess. They lead to LDA-GW ban dgap s of ZB-InN 
of 0.52 eV (Ref. EtJ) and - 0.4 eV (Ref. El obtained by 
subtracting 0.2 eV from the reported result for the W 
phase), which are smaller than the EXX result, of 0.81 
eV. 



C. Energetic position of the d bands 

our LDA results for Ed obtained using LDA-PP's, see 
Table III, are significantly higher than the correspond- 
ing experimental data, in agreement with previous LDA 
calculations. This deficiency has been attributed to the 
incomplete cancellation of the self-interaction p|. It is 
therefore tempting to speculate that the EXX calcula- 
tions, which are self-interaction free, would give improved 
d band energies. A remarkable finding of the present 
study is that this is not the case: the EXX results are 
very close to those of LDA. In GaN the EXX result is 
only 1 eV lower than the LDA results, while in InN it 
is about 0.1 eV higher. The LDA-GW approach, on the 
other hand, is found to lead to a significant downward 
shift in ener gy o f the d bands in many II- VI and group- 
Ill nitrides However, the GW results for Ed are 
systematically higher than experiment by about 1 eV. 

In order to understand the above unexpected EXX re- 
sults for Ed let us inspect the LDA and EXX eigenvalues 
of the isolated Ga, In, and N (pseudo)atoms, shown in 
Fig. 3. These results are obtained by using the same 
pseudopotentials and XC functionals as those used for 
the corresponding bulk calculations, see Sec. II. The im- 
portant features to note from Fig. 3 are as follows, (i) 
All the eigenvalues shown are shifted downward in en- 
ergy due to the removal of the unphysical self-interaction 
by the EXX approach, (ii) The downward shifts in the 
eigenvalues of the N 2s and 2p states are comparable to 
those of Ga 3d states, and larger than that of the In 
Ad states. This explains why in the GaN and InN com- 
pounds, the EXX d bands remain roughly speaking at 
the same energetic position with respect to the N 2s and 
2p derived valence bands as in the LDA. (iii) In the LDA 
calculations using EXX-PP's only the eigenvalues of the 
Ga 3d and In Ad states are appreciably shifted down- 
ward in energy, compared to the LDA calculations for all 
electrons (i.e., using LDA-PPs). It is interesting to note 
that these shifts (which are of 3.0 and 1.5 eV, respec- 
tively) closely match the downward shift of the d bands 
in the corresponding bulk calculations (see Table III). 
This leads to values for Ed in a rather good agreement 
with experiment. Similar results are obtained for some 
II- VI compounds by Moukara et al. [35|. who suggested 
that further improvements might be achieved if the EXX 
approach is applied to all electrons. The results of the 
present work show that this is not correct. 
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D. Valence band- widths 

We will focus mainly on (i) the width of the upper 
three valence bands, which mainly originate from the N 
2p states, referred to as "upper VBW". (ii) The width 
of the total valence bands (referred to as "total VBW"), 
defined as the energy difference between the top of the 
valence bands and the minimum of the lowest energy va- 
lence band which originates from the N 2s states. This 
means that the cation semicore d bands are not consid- 
ered as part of the main valence bands. 

Fig. 1 and Tables II and III show that the upper VBW 
decreases (by about 0.5 eV) by going from LDA to EXX 
calculations. The width of the lowest valence band (see 
Fig. la) also decreases, but by a much smaller amount. 
The reduction of the upper VBW is consistent with pre- 
vious EXX calculations 0, 0, . This feature can be 
explained pH as a direct consequence of the increase in 
strength of the effective potential seen by the valence 
electrons, due to the removal of the self-interaction by 
the EXX approach. The LDA-GW upper VBW's HEU 
are larger than those of EXX by about 1 eV, for both 
A1N and GaN. For a discussion about the comparison 
between the EXX and LDA-GW results see the following 
subsection. 

Fig. 1 and Tables II and III show that the total VBW is 
weakly affected by going from LDA to EXX calculations. 
The total VBW is affected through the reduction of the 
upper and lower VBW's and the difference in the down- 
ward shifts in energy of the involved bands, and these 
two contributions seem to cancel each other to a large ex- 
tent. Table II show that, unlike the case of the bandgap, 
quasiparticle corrections have sizable effects on both the 
upper and total VBW's, in the case of both LDA (i.e., 
LDA versus LDA+GW results) and EXX (EXX versus 
EXX+GW results) approaches. The LDA-GW results 
0, |4lJ for the total VBW are appreciably larger than 
those of EXX. For GaN, both the EXX and GW results 
deviate significantly from the experimental value for the 
total VBW. The reason for such bad disagreement is un- 
clear, but experimental uncertainties can't be ruled out. 



E. Further discussion 

In this subsection we discuss our EXX results in com- 
parison with those of GW, from the point of view of the 
conceptual difference between the two approaches. Fur- 
thermore, we discuss the role of using consistent pseu- 
doptentials (i.e., EXX-PP's) in the EXX calculations, 
and comment on the discrepancies between the pseu- 
dopotential and all-electron EXX results. 

We start with the comparison between the EXX and 
LDA-GW results. These two approaches provide dif- 
ferent approximations to the XC self-energy operator. 
The LDA-GW approach is a many-body technique which 
leads to approximate quasi-particle or -hole energies, 
while the EXX method is a single-particle approach. 



As shown in this work and Refs. 0, 12^) there is a 
very good agreement between the bandgaps calculated 
by both methods for wide range of semiconductors. This 
may suggest that for these systems the XC potential in 
the EXX approach (including LDA correlation) can be 
considered as a very good approximation to the LDA- 
GW self-ener gy o perator, in the optimized effective po- 
tential sense |5l|. The conduction and upper valence- 
bands states are highly delocalized states and, thus, the 
effects of charge density relaxation by the removal (ad- 
dition) of an electron from (to) these states are expected 
to be small. This is not the case for the more localized 
valence and semicore d states, as the ASCF calculations 
show 0, • This explains the good agreement between 
the EXX and LDA-GW results for the bandgaps, and 
the appreciable discrepancies between their results for Ed 
and the valence band widths. Furthermore, this also im- 
plies that one should be careful when comparing the cal- 
culated EXX (or LDA) binding energies of the semciore 
d electrons with the experimental photoemission data, 
since the latter include electronic relaxation effects. 

Now we address the question of why the EXX ap- 
proach yields highly improved bandgaps, especially for 
the conventional sp semiconductors. Following Stadele 
et al. Hj, we write the bandgap as 

E g = Ef s + A XC = Ef xx ™ + e£f + A XC . (1) 

E X X ( X ") 

Here, E g is the EXX bandgap obtained without 

any correlation, and e^f is the change in the bandgap 
due to the KS correlation potential. For conventional sp 
semiconductors it has been shown that E g rs E g XX ^ , 
and t^f. is quite small at the LDA or GGA levels (~ 0.1 
eV) . These results have led Stadele et al. [21| to conclude 
that Axe is quite small, although the corresponding dis- 
continuity in the EXX potential alone is very large, com- 
pared to the bandgaps of these systems |2l| . On the other 
hand, the Random-phase approximation (RPA) investi- 
gations [lli|53|,|54j have shown that the KS bandgaps de- 
rived from these GW calculations are very close to those 
of LDA, for bulk Si and C. This means that RPA is 
quite large. Thus, according to these RPA calculations, 
the EXX approach provide improved bandgaps because 
of a large cancellation between Axe and the errors in 
the LDA (or GGA) e£f . 

The effects of employing consistent pseudopotentials in 
the EXX calculations (on the calculated bandgaps and 
valence band-widths) can be easily seen by comparing 
the EXX or LDA results obtained using both LDA- and 
EXX-PPs, listed in Tables II and III. The interesting fea- 
tures to note are: (i) if the same set of pseudopotentials 
(LDA or EXX) is used in both calculations, the band- 
structure modifications by going from LDA to EXX cal- 
culations are almost independent on the employed pseu- 
dopotentials. (ii) For the three nitrides considered, these 
effects are quite small on the above properties, contrary 
to Ed, expect for the total VBW of InN obtained with 
the semicore d electrons included as core states. This 
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can be easily understood from the shifts in the eigen- 
values of the pseudoatoms shown in Fig. 3. It should 
be noted that this is not always the case: the use of 
EXX-PPs in the LDA calculations leads to an increase 
in the bandgap of GaAs by about 0.6 eV, and to a cor- 
rect ordering of the conduction band minima of Ge |2l| . 
(iii) However, the advantages of using consistent pseu- 
dopotentials in the EXX calculations of the considered 
nitrides are clear: the comparison between EXX results 
obtained using EXX-PPs with experiment is better than 
those obtained employing LDA-PPs, especially for the 
the bandgap of GaN and Ed of both GaN and InN. 

After the completion of this work, all-electron full- 
potential EXX results have been reported jEH- These 
results are considerably different from those of pseudopo- 
tential EXX calculations, especially for the bandgaps of 
sp semiconductors and the energetic position of the d 
bands. For the latter the all-electron EXX results are in 
very good agreement with experiment, contrary to our 
results (see above and Ref. |24J). These discrepancies 
have been attributed mainly to the core-valence interac- 
tion |5^. We strongly believe that such a conclusion is 
misleading because of two reasons: (i) We have shown 
(in Sec. III.B) that our pseudopotential EXX results for 
Ed can be easily explained in terms of the eigenvalues 
of the corresponding isolated atoms, where core- valence 
interaction is not an issue, (ii) There are significant dif- 
ferences in the implementation of the EXX scheme in the 
two approaches (see Refs. |^,|5^for details), which may 
account for the above discrepancies. However, the reason 
behind these discrepancies can be clarified by performing 
pseudopotential EXX calculations with the entire semi- 
core shell included as valence, which is currently under 
consideration. 



IV. CONCLUSIONS 

The exact-exchange (EXX) Kohn-Sham density- 
functional theory is used to calculate the electronic struc- 



ture of the zinc-blende (ZB) form of A1N, GaN and InN, 
with and without including the cation semicore states. 
The effects of using consistent (i.e., EXX) pseudopoten- 
tials in the EXX calculations are also investigated. The 
changes in the bandgap and energetic position of the shal- 
low d bands by going from the ZB to wurtzite (W) forms 
have been studied, at the LDA level. Our EXX results 
for A1N and GaN, obtained with the semicore d elec- 
trons treated as part of the core, are found to be in good 
agreement with similar previous EXX calculations, GW 
results and experiment. Treating the semicore d electrons 
as valence leads to a large reduction in the calculated 
EXX bandgaps, and the results are in excellent agree- 
ment with the pseudopotential (when the entire semi- 
core shell is included as valence) and all-electrons GW 
results. This may indicate that the inclusion of the en- 
tire semicore shell as valence is not necessary in the EXX 
and EXX+GW approaches. Contrary to common belief, 
the removal of the self-interaction by the EXX approach, 
does not account for the large discrepancy for the posi- 
tion of the semicore d bands between the LDA results 
and experiment. Most of the opening of the bandgaps 
of the considered nitrides in the EXX approach is due to 
its application to the valence electrons, and its applica- 
tion to the core electrons has pronounced effects on the 
position of the d bands. The EXX method gives an ap- 
preciable bandgap for InN when the semicore d electrons 
are included as valence, which provides the basis for more 
accurate GW calculations. 
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Table II: Selected EXX and LDA electronic structure 
properties of the ZB phase of the nitrides considered, 
obtained with the semicore d electrons treated as core 
states. All tabulated quantities are in eV. 





LDA 


EXX 


LDA-GW EXX-GW 


A1N 


4 
y 


4.27 a , 


5.99° 


Q QC,d 


5.98 c 




4.03 b 


5.76 b 








4 20 c,rf,e 


5.66 c ' e 








3.27 a 


5.18 Q 


49 c,d 


4.89 c 




3.16 b 


5.08 b 








3.20 c ' d 


5.03 c ' e 








3.24 e 








upper VBW 


5.91 a 


5.27 a 


g 7 c,d 


6.29 c 




5.96 b 


5.33 b 








60 cd 


5.24 c 






total VBW 


14.81° 


14.59° 


17.00 c - d 


15.32 c 




15.16 b 


14.95 b 








15.10 c ' d 


14.85 c 








14.74 6 


14.86 6 






GaN 


4 


2.20 a 


3.68 a 


3.59' 


3.44 c 




1.97 b 


3.52 b 


3.10 c - d 






2.00 c 


3.49 c ' e 


2.76 9 






2.10 d 










1.9 e 










1.88 s 








upper VBW 


6.85° 


6.30 Q 


6.88* 


7.17 c 




6.82 b 


6.23 b 


7.6 C 






6.8 C 


6.23 c 


7.8 d 






7.4 rf 








total VBW 


15.85 a 


15.35 a 


16.70 c 


16.05 c 




15.82 b 


15.63 b 


17.80 d 






15.50 c 


15.64 c ' e 








16.30 d 










15.75 6 








InN 




0.17 a 


1.89 a 


1.31" 






-0.13 b 


1.49 b 






upper VBW 


5.52 Q 


4.88 a 








5.67 b 


5.09 b 






total VBW 


13.98 a 


13.82° 








14.50 b 


14.37 b 







"Present work, using LDA pseudopotentials. 
b Present work, using EXX pseudopotentials. 
c Ref. [8]; d Ref. [41]; e Ref. [21]; ^Ref. [10]; 
fRef. [42] ; h Ref. [27]. 



Table III: Selected EXX and LDA electronic structure 
properties of the ZB phase of GaN and InN, with the 
semicore d electrons treated as valence states. For com- 
parison some results for the W phase are shown, denoted 
by the superscript w. All tabulated quantities are in eV. 





LDA 


EXX 


LDA-GW 


Expt. 


GaN 




1.65° 


2.67 a 


2.88 c 


3.2 d 




1.66 b 


2.88 b 


3.03™' 9 


3.3 e 


E d 


-13.8 Q 


-12.2 a 


-15. T 


-17.7 h 




-16.9 b 


-14.8 b 


-16.4^'f 






-13.8<"^ 








upper VBW 


7.32° 
7.03 b 


7.03 a 
6.63 b 


7.33 c 




total VBW 


16.00 a 
16.04 b 


16.04 a 
16.00 b 




14.2 d 


InN 




-0.43 a 


0.61° 


~ 0.0™ J ' 9 


0.7-1.0 , " J 




-0.39 b 
-0.36* 


0.81 6 




l.9 w ' k 


E d 


-13. 5 Q ' 1 


-12.2 a 




-U.9 W ' 1 




-14.9 b 


-13.4 b 




-16.7 w ' m 












upper VBW 


6.10 Q 
5.85 b 


5.82 a 
5.48 b 






total VBW 


14.60 a 
14.71 b 


14.67 a 
14.73 b 







"Present work, using 
b Present work, using 
c Ref. [lfj]; d Ref. [45J; 
h Rcl [48]; l Ref. [27]; 
m Ref. [5J. 



LDA pseudopotentials. 
EXX pseudopotentials. 

e Ref. [46J; f Ref. [4JJ; s Ref. [11]. 

J'Ref. [28]; fc Ref. [29]; z Ref. [49]. 




FIG. 1: Band structure of GaN calculated without (a) and 
with (b) treating the Ga 3d electrons as valence states. The 
EXX pseudopotentials have been used in both LDA and EXX 
calculations. 
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FIG. 2: The atomic EXX and LDA all-electrons potentials 
of Ga. The right asymptotic — 1/r behavior of the EXX po- 
tential should be noted. The EXX and LDA eigenvalues of 
the Ga 3s and 3p are shown. 
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FIG. 3: Eigenvalues of the Ga, In and N (pseudo)atoms. 
Dashed lines: LDA calculations using LDA-PP's. Dotted 
lines: LDA calculations using EXX-PP's. Solid lines: EXX 
calculations using EXX-PP's. 



